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ABSTRACT 

We present a circular cross-correlation tests for the phases of the Internal Linear Combination Map 
(ILC) and WMAP^s foregrounds for all K-W frequency bands at the range of multipoles £ < 50. We 
have found significant deviations from the expected Poissonian statistics for the ILC and the foregrounds 
phases. Our analysis shows that the low multipole range of the ILC power spectrum contains some of 
the foregrounds residues. 

Subject headings: cosmology: cosmic microwave background — cosmology: observations — methods: 
data analysis 



related with foreground contaminations (see also Chiang 
et al. 2003) produced most likely by the Galactic dust 
emission. 



2. ILC AND FOREGROUNDS PHASE CORRELATIONS 

The basic idea of the phase analysis of the CMB maps 
was proposed and discussed in Chiang and Coles (2000), 
Naselsky, Novikov and Silk (2002), Chiang et al. (2002), 
Chiang et al. (2003), Naselsky et al. (2003), Matsubara 
(2003) 

The CMB signal from can be expressed as a sum over 
spherical harmonics: 

oo I 

aim = \aim\ exp(i^'£,„), 



(1) 



1. INTRODUCTION 

The recently-published Wilkinson Microwave 
Anisotropy Probe ( WMAP) data sets (see Bennett et 
al. 2003 a-c, Hinshaw et al. 2003 a-b)) strongly promote 
the development of the high sensitive statistics for testing 
of the properties of the derived foregrounds and fore- 
ground cleaned maps (Komatsu et al. (2003), Tegmark, 
de Oliveira-Costa and Hamilton (2003), Chiang et al. 
(2003), Dineen and Coles (2003), Park (2003), Gaztanaga 
et al. (2003), CoUey and Gott (2003), Naselsky et al. 
(2003)). The WMAP team produces the Internal Lin- 
ear Combination (ILC) map and maps for the foregrounds 
(synchrotron, free- free, dust emission) for each K~W bands 
which are the basis for our analysis. Homogeneous and 
isotropic CMB Gaussian random fields, as a result of 
the simplest inflation paradigm, is a crucial test for the 
early stages of the cosmological expansion. Because of the 
pronounced non-Gaussianity of the foregrounds it would 
be natural to expect that possible non-Gaussianity of the 
ILC and other cleaned maps derived from the WMAP data 
would reflect directly contamination of the foregrounds at 
different levels. 

In this Letter we present the result of a statistical test of 
the coupling the ILC map and the foreground taking from 
the WMAP data^. The test is based on the circular cross- 
correlations analysis of the ILC and foregrounds phases. 
Our method complementing the mentioned above ones, 
exploits a natural assumption that phases of the "true" 
CMB signal should not correlate with the phases of the 
foregrounds^. This allows us to detect significant ILC-"W 
band foregrounds" cross-correlations. These peculiarities 
determine the accuracy of the power spectrum estimation 
for the ILC map and the statistical properties of the ILC 
signal. We would like to point out that the detected non- 
Gaussianity of the ILC map (see for example, Eriksen et al. 
2003, Vielva et al. 2003, Coles et al. 2003) is most likely 
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where aim are the coefficients of expansion, \aim\ is the 
modulus and "^im is the phase of each aim harmonic. 
Homogeneous and isotropic CMB Gaussian random fields 
(GRFs), as a result of the simplest inflation paradigm, 
possess Fourier modes whose real and imaginary parts are 
independently distributed. The statistical properties of 
GRF are completely specified by its angular power spec- 
trum Ci, 



C{e,e\m,m') = {aima*i,m') 

(|a£m||a^'m'|)(exp (i(*i?™ - '^tm'))) = C'(^)^££'^mr 



(2) 



where (..) means averaging over the modulus and phase 
distribution functions. In Fig.l we plot the cross- 
correlation of phases between the ILC map and the com- 
mon foregrounds defined as a sum over synchrotron, free- 
free and dust emission maps for each frequency band^. 
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Flc;. 1. — The cross-correlation of phases for ILC and the K-W 
bands foregrounds. All the foreground maps and derived phases 
were taJjing from the WMAP web-site. 

As one can see from Fig.l the distribution of the ILC 
map phases and corresponding foregrounds seems to be 
non-correlated one. However, in order to test such a dis- 
tribution, for these phase diagram we will apply circular 
statistics (Fisher 1993) to investigate the possible cross- 
correlation quantitatively. 

3. CORRELATIONS OF THE ILC AND FOREGROUNDS 
PHASES 

Nasclsky ct al. (2003) have shown that for linear meth- 
ods of the foreground separation from the CMB signal we 
can expect to find some cross correlations between the 
phases of the cleaned signal, and the foregrounds, V/, 
which can be used to characterize the degree of separa- 
tion achieved. To do this, here wo will use the circular 
correlation coefficients which reflect directly the 27r-phase 
periodicity. 

For the K-W frequency band the maps were taking from 
the WMAP web site* and all the phases are obtained by 
the spherical harmonic decomposition using the HEALPix 
(Gorski et al. 1999) and the GLESP (Doroshkevich et 
al. 2003) codes. We consider the ranges of multipoles, 
2 < £ < £max = 50 for which the WMAP own foreground 
are presented in the WMAP web site. In addition to the 
WMAP foregrounds for each band we produce the five 
maps as a difference between signal {S) in the band and 
the ILC signal: F = S - ILC. We will call the F-map 
as the derived foreground. Let "tp and 4> be the foreground 
and ILC phases, respectively, for given value of £ and all 
the corresponding values of m. Following Fisher (1993) we 
define the statistics 

Xm = COs(V'm), Um = Sin(V'm), 

l^m = cos(0m), Um = sin((/)m), 
* http: //lambda. gsfc .nasa.gov/product/map/m_products . cfm 



Mp = {l/£ 

max 

m 



((/)) = tan~^( 



MU 



) XI exp {ip{il^, 



(V) = tan-i( 



(4 



(3) 



where Mp and M\p are the p-th trigonometric moments 
of the phase samples, (0) and {tp) are the correspond- 
ing mean directions of the samples, Rsf{^) is the circular 
cross-correlation coefficient in each mode £ and r^/ is the 
mean circular cross-correlation coefficicmt for all phases. 
For m = and for all £ phases (f){£,0) = ^^{£,0) = and 
here we neglect them. 

For the K-W bands the circular coefficients, Rsf{() are 
plotted in Fig. 2 for the ILC and its own foreground, and 
for the ILC and the derived foreground. 
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Fig. 2. — The circular correlation coefficient between the cleaned 
signal and the foreground phases in V— W cliainicls vs. the harmonic 
index, I. Thick solid line corresponds to the WMAP own fore- 
ground. Thin solid line corresponds to derived foreground. Shadow 
area represent la error bars level taking from 200 random realiza- 
tions. 

As is seen from the Fig. 2, for the first three channels 
these coefficients are quite moderate and do not exceed 
the random scatter (Icr) obtained from 200 random real- 
izations. Note that, for all the bands, the shape of the 
functions Rsf{^) are quite similar to each other which re- 
flects the strong correlation of phases in all the channels 
(Naselsky et al. 2003). As one can see from Fig. 2, the 
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cross-correlation of the ILC phases and the derived fore- 
ground seems to be more significant. The same tendency 
fohows from the estimation of the mean coefficients, rgf, 
for 2 < < 50, listed in Table 1 for the two samples of the 
foregrounds. As is expected, for the first three channels 
these coefficients are small but they become significant for 
channels V and W. 

Table 1 

Circular cross-correlation coefficients between phases of ILC 
and the WMAP own foregrounds (ILC*°') and the ILC and 
derived foregrounds (ILC''*'). 





K 


KA 


Q 


V 


W 




-0.026 
-0.017 


-0.031 
0.002 


-0.030 
0.022 


-0.033 
0.112 


-0.033 
0.262 



To test the sensitivity of circular statistics, we used the 
so called "WMAP simulator" , in which we simulated the 
W band with all own WMAP foregrounds and the CMB 
signal with random phases instead of the ILC signal^. Af- 
ter that we repeat the circiilar phase analysis in the same 
way as it was done for the ILC map. Our results are shown 
in Fig. 3. 
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Fig. 3. — Raf for the "WMAP simulator". Thick solid line corre- 
sponds to the simulated random CMB signal, cross-correlated with 
WMAP own foreground. For the CMB and derived foregrounds 
we have the same curve. The dashed line is the cross-correlation 
between the foreground, defined as the W band signal minus the 
simulated CMB signal (not the ILC !), and the simulated CMB 
signal. 



For the simulator, corresponding mean cross-correlation 
coefficients are r^/ = —0.036 for the CMB-Foreground 
correlations and Tsf = —0.194 for "wrong" foreground 
(actual signal minus simulated CMB), corresponds to the 
dashed line in Fig. 3. Thus, if the CMB map and fore- 
ground are extracted correctly, corresponding mean cross- 
correlation coefficient should be less then statistical error 
(~ 5%), while the wrong foreground and CMB should have 
significant positive or negative bias. For the ILC map it 
means that a positive bias for fg/ for the V and W bands 
(see Table 1) reflects directly the cross-correlation with the 
foreground phases. 

4. CONCLUSION 

The circular cross-correlation analysis for the ILC map 
and K-W foregrounds and shows that the ILC map has sig- 
nificant correlation with the derived foregrounds obtained 
by subtraction the ILC map from the V and W bands 
signals. Some of this correlations could be linked with 
non-Gaussianity of the ILC map detected by Park (2003), 
Eriksen et al. (2003), Coles et al. (2003), Vielva et al. 
(2003). For example, Vielva et al. (2003) has reported 
about characteristic scale of non-Gaussianity in order to 
10°, which corresponds to ^ ~ 10. Our analysis shows 
that the cross-correlation coefficients Vsf at ^ = 11 have 
the highest maxima for both V and W bands. At high mul- 
tipole range £ ~ 40 corresponding peculiarities are above 
95% confidential level. If wc take into account that for 
£ ~ 40 all m ^ 35 — 40 lies along Galactic plane in the 
map, we can conclude that they are related with the V 
and W bands signal at the same longitude. Some of the 
peculiarities need an additional investigation which would 
be published soon. 
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